Absence of domain wall roughening in a transverse-field Ising model with long-range 

interactions 
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We investigate roughening transitions in the context of transverse-field Ising models. As a modifi- 
cation of the transverse Ising model with short range interactions, which has been shown to exhibit 
domain wall roughening, we have looked into the possibility of a roughening transition for the case 
of long-range interactions, since such a system is physically realized in the insulator LiHoF4. The 
combination of strong Ising anisotropy and long-range forces lead naturally to the formation of 
domain walls but we find that the long-range forces destroy the roughening transition. 

PACS numbers: 75.30.Kz, 75.60.-d 
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I. 



INTRODUCTION 



Magnetic systems present us with the opportunity to 
study not only classical but also quantum critical phe- 
nomena and thus provide us with unique insights in con- 
densed matter physics and the intertwining of classical 
and quantum mechanics^. Of particular interest is the 
transverse field Ising modeli*^, whose representation in 
terms of Pauli spin matrices is 



N 



N 



(i) 



where h x represents an applied transverse magnetic field 
and Jij are coupling constants. Notice that in the ab- 
sence of the magnetic field the Hamiltonian is diagonal 
in the a z basis, and the system is simply a classical Ising 
model. As the magnetic field is turned on, a x opera- 
tors are introduced which do not commute with the a z 
operators. Thus, turning the transverse field on or off 
essentially turns quantum mechanics on or off in the sys- 
tem. The spins of this system align at temperature T=0 
in a ferromagnetic ground statei whereas at high temper- 
atures the system becomes disordered. In the absence of 
a transverse magnetic field, a continuous phase transition 
occurs between the paramagnetic and ordered ferromag- 
netic states. This phase transition is driven by thermal 
fluctuations. 

Applying a transverse magnetic field, h x , perpendicu- 
lar to the axis of preferred magnetization, can also cause 
a transition between the ferromagnetic and the disor- 
dered states even at zero temperature^. This behav- 
ior is driven by quantum zero-point fluctuations of the z 
component of the spins due to the transverse field. 

We now imagine the system described by (1) has do- 
main walls - this is achieved in principle by imposing 
appropriate anti-periodic boundary conditions. It is ex- 
pected that, in the absence of the transverse field, at 
T=0 the domain walls would be flat. As we increase 
temperature we would first observe nucleation of steps in 
the interface. At some temperature, called the roughen- 
ing temperature Tr, the entropy of thermal fluctuations 




FIG. 1: LiHoF4 behaves essentially as an Ising model with 
long range interactions. Interesting physics arises with the 
application of a transverse magnetic field. 



of the interface would dominate the interfacial energy 
and the interface would become rougbA At tempera- 
tures above the roughening transition, the amplitude of 
fluctuations in the surface's position scales as L e where L 
is the system size and is an exponent characterizing the 
nature of the interface. Thus the interface fluctuations 
would diverge as the system length increases to infinity. 
However, we note that if 6 < 1 then the system size would 
diverge faster than the fluctuations and so the fluctua- 
tions of the interface would not overwhelm the bulk order 
in the system. This would allow for a roughening transi- 
tion to take place at a temperature lower than the bulk 
order-disorder phase transition temperature. 

A physical realization of the transverse Ising model 
as described by Eq. Q is provided by the insulating 
magnet 2,3 LiHoF4. In this case, the coupling constants 
Jij are actually no longer near-neighbor but dipolar in 
nature and so decay as the inverse distance cubed. 
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In LiHoF4 the Ho 3+ ions are responsible for the magnetic 
behavior. The crystal field splitting of the Ho 3+ states 
is such that the ground state is doubly degenerate and 
well below the higher states, thus leading to a very strong 
Ising anisotropy. The two spin states for each Ho 3+ ion 
point along the z-axis, as shown in Fig. ^ The phase 
diagram, Fig. [21 is indicative of a quantum phase transi- 
tion at T = and at a critical transverse magnetic fields 
H x = H c ^49 kOe, as well as a classical phase transi- 
tion between ferromagnetic and disordered states in the 
absence of a magnetic field at T c =1.53 K. 

The dipolar long-range forces and the insulator nature 
of LiHoF4 lead to the natural formation of needle-like al- 
ternating domains of antiparallel magnetization, without 
the need to force antiperiodic boundary conditions upon 
an experimental sample. Since LiHoF4 is an insulator, 
the electrons and hence the constituent spins of the sys- 
tem are localized. The large anisotropy of the system 
means the domain walls which separate regions of oppo- 
site magnetization in a LiHoF4 sample are sharp and well 
defined, in the sense that as we move across a domain in- 
terface the system abruptly switches from one magnetiza- 
tion state to another. This is to be contrasted to domain 
walls found in metals which are extended and provide for 
a continuous smoother transition of the magnetization as 
we move from one domain to the next. Thus, the first 
impression we get is that LiHoF4 appears as an ideal sys- 
tem to study roughening transitions of domain interfaces. 
However, our analysis suggests otherwise. 

In what folllows we first briefly review the case of only 
short-ranged exchange interactions being present, where 
we would expect to see a roughening transition, as de- 
scribed by Fradkin 5 . Then, we demonstrate that in the 
case of dipolar forces such as seen in LiHoF4 the long- 
range nature of the interactions drives the roughening 
temperature up to the bulk transition temperature, so 
that the domain walls remain flat throughout the ferro- 
magnetic regime. We show how to obtain a Hamiltonian 
for a field theoretical description of a single domain wall, 
which results in a modified sine-Gordon model that be- 
haves as the regular sine-Gordon model but with effec- 
tively higher dimensionality. We see how this behavior 
is verified by a renormalization group calculation using 
a smooth cutoff and conclude that the effective higher 
dimensionality of the domain wall, d>2, indicates the 
absence of a roughening transition and the persistence 
of flat domain walls in the ferromagnetic phase. Finally, 
we mention some further considerations for stepped in- 
terfaces. 



II. TRANSVERSE ISING MODEL WITH 
ROUGHENING: SHORT RANGE 
INTERACTIONS 

Using an Ising model for a two-state spin system with 
appropriate boundary conditions, a two-region problem 
can be set up with an interface of finite width separating 
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FIG. 2: LiHoF4 undergoes both a classical thermal phase 
transition as well as a quantum phase transition with the 
application of a transverse magnetic field, after Ref. |2|. 

regions of opposite spin. Investigations of this interface 
for the case of short-range interactions show that it un- 
dergoes the roughening transition&iSuLiSiil. A numer- 
ical estimate for the roughening transition was carried 
out some time ago by Weeks et al. 12 for a 2-d interface 
in a classical Ising model, obtaining a roughening criti- 
cal temperature of 0.57T C , where T c is the bulk critical 
temperature for the order-disorder transition. It would 
be useful to redo this calculation with higher precision 
using newer computational techniques. 

In the case of two-dimensional interfaces it was pro- 
posed by Fisher and Weeks^ 3 - that the interfaces of three- 
dimensional quantum crystals are always smooth at zero 
temperature. This suggests that there can be no quan- 
tum roughening transition. Fradkin 5 * investigated the 
suggestion further and verified it for two different models, 
namely, a model describing the solid to vacuum interface 
of a three-dimensional quantum crystal, which is similar 
to the model used by Fisher and Weeks, and a model 
describing the interface of the three-dimensional quan- 
tum transverse Ising model. Fradkin obtains the effective 
Hamiltonian of the two-dimensional interfaces, both for 
the quantum crystal model, 

Hqc = -A'^cos[p(r) -p(r')] 

r 

+ ~ E Hr)-n{r')]\ (2) 

<r,r '> 

and for the transverse Ising model, 
Hti = — cosp(r) 

r 

+ \ E H^)-n{r')]\ (3) 

<r,r ' > 

where p is canonically conjugate to n, the discrete height 
variable. 
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Notice here the difference in the kinetic energy term 
which arises from the difference in symmetry between 
the different models and gives rise to different dynamics 5 . 
For the quantum crystal model the particle number N — 
^-n(r) is conserved while for the analogous quantity for 
the transverse Ising model it is not. Since our model for 
LiHoF4 is the transverse Ising model we take into account 
the considerations by Fradkir. 5 instead of the Fisher and 
Weeks modelii 

The 2+1 dimensional quantum interface model for the 
transverse Ising is dual to a three-dimensional general 
Coulomb gas which is known to have only a conducting 
plasma phasoii Since the surface has to be smooth at 
T=0 this implies that roughening can only occur if the 
bulk loses its long-range order, so the critical magnetic 
field parameter for the roughening and order-disorder 
transitions will coincide, i.e. Hr = h c . 

For finite temperatures, where the transition is classi- 
cal, we can ignore the cosine kinetic energy terms since 
these become unimportant, and notice that the two mod- 
els in Eqs. J2J and J2J become the same, namely a realiza- 
tion of the discrete Gaussian SOS model 5 . A continuum 
version of the discrete Gaussian SOS model is the sine- 
Gordon model 15 . Briefly, the sine-Gordon model hamil- 
tonian contains two terms: 



2-6 Interface 



H sG = I d 2 r 



27m (r) 



- (Vn(r)) 2 - ycos 
2 \ a 



(4) 



The first term describes the cost of local fluctuations of 
the interface. The second term takes the effects of the lat- 
tice into account, and the fact that the interface height 
field n(r) can really only take discrete values, namely 
multiples of the lattice spacing a. This periodic term pro- 
motes the pinning of the interface onto equidistant par- 
allel planes, and thus encourages the interface to remain 
smooth. The sine-Gordon model undergoes a Kosterlitz- 
Thouless phase transitio n 15 : 16 and displays two distinct 
phases corresponding to the smooth surface at low tem- 
perature, where the cosine term is perturbatively impor- 
tant, and a rough phase at high temperatures, where the 
cosine term is unimportant. Thus we expect a rough- 
ening phase transition at some finite temperature, Tr, 
less than the bulk order-disorder transition temperature. 
Given that the interface has to be smooth at T = 
we propose the phase diagram shown in Fig. [21 where 
we include, for completeness, the related renormalization 
group (RG) flows. For the two-dimensional interface we 
notice that in the temperature region Tji < T < T c the 
flows are opposite for the two different phase transitions. 
This is related to the fact that the onset of high fluctua- 
tion for the interface occurs at a temperature below the 
critical order-disorder transition where the bulk system 
loses its long-range order. 

It is also interesting to consider a one-dimensional in- 
terface in the transverse field Ising model, in which case 
we expect the domain wall to be rough at any finite tem- 
peratures. However, as pointed out by Fradkin 5 -, at zero 
temperature the transverse Ising model is dual to a two- 
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FIG. 3: For a two-dimensional interface we expect a thermal 
roughening transition at some finite temperature. At zero 
temperature the critical transverse field for the roughening 
transition coincides with that of the bulk order-disorder tran- 
sition. The two sets of arrows indicate the expected qualita- 
tive renormalization group flows for the different phase tran- 
sitions: the bulk order-disorder transition (plain arrows) and 
the interface roughening transition (double-headed arrows). 



dimensional Coulomb gas 17 . This is dual in the contin- 
uum limit to a two-dimensional sine-Gordon model and is 
known to have a metal- insulator phase transition 14 ' 16,18 
at some value of h x at T=0: we expect a smooth inter- 
face for small h x ; as we increase h x we go through a phase 
transition to a rough interface for some tiR, before the on- 
set of the order-disorder transition at h c . The proposed 
phase diagram and qualitative RG flows are shown in 
Fig. 0] We see that for the roughening transition, tem- 
perature is always a relevant parameter and flows to a 
fixed point at infinity. 



III. LiHoF 4 TWO-DOMAIN SYSTEM: 
LONG-RANGE INTERACTIONS 

As mentioned in the introduction, the behavior of 
LiHoF4 can be understood in terms of a quantum Ising 
model in the presence of a transverse magnetic field^, 
in which the states of the Ho 3+ ions are represented by 
the Ising spins | |) and | J,). The effective Hamiltonian 
looks likei 



H = 



I 

2 



(5) 



<i,3> 



where h x corresponds to an applied transverse magnetic 
field and J ex is an effective antiferromagnetic exchange 
term that is added to obtain an effective Hamiltonian 
that matches the experimental results. 
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FIG. 4: For a one-dimensional interface the system is rough 
at all finite temperature. The roughening transition oc- 
curs at zero temperature. In addition we expect a quantum 
roughening phase transition. The two sets of arrows indi- 
cate the expected qualitative renormalization group flows for 
the different phase transitions: the bulk order-disorder tran- 
sition (plain arrows) and the interface roughening transition 
(double-headed arrows) . 



Since LiHoF4 forms domains naturally, we consider for 
our calculations a system consisting of two semi-infinite 
three-dimensional domains of antiparallel magnetizations 
(directed in the positive or negative z-directions) that 



are separated by a domain wall set at x = 0. We want 
to calculate the configurational energy for fluctuations 
in the domain wall with respect to a flat interface. We 
proceed by assuming a height for the domain interface to 
be a general function ip(y,z), so that the magnetization 
vector for the system becomes, 



A? = {0, 0, m sgn(x - ip(y, z)}. 



(6) 



We consider the classical magnetostatic problem and af- 
ter we relate a "magnetic charge density" q{x, y, z) to the 
magnetization, 



q{x,y,z) 



-V-M, 



(7) 



we see that in analogy with electrostatics, the extra en- 
ergy associated with having a non-flat profile is given in 
the continuum limit byiS 



(8) 



where <f> is the magnetic potential given in three dimen- 
sions by 



q(f>) 



(9) 



This results in the following expression for our system 



U = 2m, 



/ dydzdy'dz' 



d z ^(y,z)d z >ip(y\z') 



My, z) - ip(yi, z')Y + (y- y') 2 + (z- z'Y 
I 



1/2- 



(10) 



We now assume that the profile of the interface 
does not vary greatly from one position to another 
and so we can expand the denominator to zeroth or- 
der in [ip(y, z) — %p(y' , z')], which is certainly true in the 
"smoother" regime of 9 <C 1 and probably asymptotically 
correct for all 6 < 1. We can rewrite the energy in this 
approximation in Fourier space to obtain: 



Un = Airmi 



dkydk 2, 



{2-Kf \k\ 



\?P(k v ,k z )\ : 



(11) 



Odd order terms vanish and higher even-order terms do 
not display the singular behavior of the zeroth-ordered 
term which will ultimately dictate the critical behavior 
of our system. 

A point of concern here is that, in the continuum model 
used, the fluctuations in the y-direction do not have any 
energy cost since they do not change the angle between 
magnetization up and down in the z-direction as we move 
across the interface. Thus we should consider fluctua- 



tions in the y-direction in a lattice model rather than in 
the continuum to see their effect. A further complication 
arises, since in addition to the dipolar Ising model inter- 
action which would raise the energy of the system there 
could be an effective exchange interaction in LiHoF4 as 
indicated in Eq. (J5J, tending to lower the energy as the 
domain wall fluctuates. This effective exchange interac- 
tion has recently been discussed by Chakraborty et a£, 
who obtained an estimate of J ex from Monte Carlo sim- 
ulations, 



0.75 J = 0.75 x 0.214K". 



(12) 



To get a feel for the system stability to deformations 
and the domain wall surface tension we took a LiHoF4 
system starting with two semi infinite domains of anti- 
parallel magnetization (a cube with sides extending to 
± N) and then we introduced a semi-infinite unit step 
deformation in the chosen y-direction. Briefly, this is 
a lattice sum calculation where the dipolar interactions 



5 



were summed over the infinite Bravais lattice, which 
for LiHoF4 has tetragonal unit cells with dimensions 
(a, a, c) = (1, 1, 2.077)a, where a = 5.175A is the lat- 
tice constant. Each unit cell has four spins, at the posi- 
tions {0,0,0},{0,f, §},{§,§, -§} and {|,0,-|}. The 
sums show slow convergence, ~ \, so we used Ewald 20 
summation to obtain faster convergence. Our results in- 
dicate that the deformation causes an increase in en- 
ergy from the dipolar term of \ x 2(27V + 1) x 0.927J 
whereas the exchange interaction would cause a decrease 
by ix2(27V+l)x0.75J. Thus, we have an overall positive 
energy cost to a single step deformation and a positive 
surface tension associated with this deformation. Taking 
the zeroth order term of the Hamiltonian Eq. I jlljl and 
following the above discussion of our numerical consider- 
ations for the deformation of an interface in equilibrium, 
we add an extra tension term, and consider the free in- 
terface to be modelled by 



U' Q = 4nml 



d k ( k ~ 9 



(13) 



We note that at this Gaussian level the individual modes 
are not coupled and hence we can evaluate the average 
of Uo assigning a variance of m to each mode. Thus we 
can write down the mean correlated height difference, G, 
as 



G = ((^(y,z)-i/j(0,0)) 2 ) 



d 2 k 



l7r ) 2 4irf3m 2 (f| + jak 2 



(14) 
(15) 



In the large distance limit, G ~ O(L ), thus indicating 
that the interface has bounded fluctuations and is smooth 
at all temperatures. 



A. Long-range sine-Gordon model 

We now wish to investigate any possibility of critical 
behavior, by performing a renormalization group analy- 
sis. Starting from the zeroth order term in our modified 
energy expression Eq. (|13fl we now include a sine-Gordon 
term to take into account the existence of the lattice and 
write the full action as 



(3H V 



d 2 k 



\k\ 



-yak 2 ip(k)i/)(-k) 



Qu / ,2 

— / a rcos 



V a 



}■ 



(16) 



where (3 = fcgT, a is a lattice constant, and 7 and g u 
are couplings for the surface tension term and fugacity 
of the system respectively. To make the coupling con- 
stants dimensionless we take pmg — > ?\,f3g u — > g, and 
in addition we define a dimensionless field in real space, 



<p(k) = ^tp(k) so that our action becomes 



+ 



d 2 k 
(2tt) 2 



k 2 

— +-yk 2 
a\k\ +7 



d 2 r cos(<j>(r)). 



(17) 



The above action is just the familiar sine-Gordon field, 
with an extra coefficient k 2 /a\k\ in the quadratic part, 
which accounts for the long range nature of the inter- 
actions. The sine-Gordon model can be renormalized 
via a smooth cutoff approach 1 - . Note here that if we 
now try to renormalize the above action using standard 
isotropic rescalings then the end result is that the singu- 



lar - 



a\k\ 



term cannot be renormalized and will dominate 



- since the renormalization of non-singular couplings can- 
not generate singular terms. Therefore, let us consider a 
different rescaling for the x and y directions while, as is 
usual, also requiring that the field <fi(r) does not change 
under rescaling so as to preserve the lattice structure. We 



rescale so that z' 



— bs 

e z, y 



s y, K 



e as k y , cf)'(k r ) = e-( b+a >(j){k), and g' = e {b+a >g, where 
a,b > 0. 

We concentrate on the terms quadratic in the field and 
choose a and b so that the two terms in the quadratic 
prefactor scale the same way under the transformation. 
This cannot be done consistently unless we assume a < b 
which leads to the choice a = |6. We choose a = 1 and 
notice that terms which behave as k 2 become irrelevant 
in this RG scheme. 

Bearing in mind the above discussion and in addi- 
tion introducing a smooth momentum cutoff function 
f(k y ,k z ) in the simplified action to take the effects of 
the lattice into account, we finally obtain the rescaled 
action 
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d 2 k'ei 



,k z )\a\k' y \ 



IT J (27T) 2 f(k y , 

+l(k l2 )}W{k')\ 2 + ¥L J d 2 rcos(0(r))(18) 

This has the form of the original hamiltonian if we define 
our new couplings as: 

9' = e^g, (19) 
m' 2 = e^m 2 , (20) 

and the new cutoff function: 

k' z ) - f(k y , k z ) = f(e~ s k' y ,e-^k' z ). (21) 

We try to pick the cutoff function in a way consistent 
with the spatial anisotropy of the problem, and also for 
future convenience. One reasonable choice is 



f(k y ,k z ) =e -« 2 (^+7a|fe,| 3 ). 



(22) 



The RG calculation (described in more detail in the 
Appendix) now proceeds as usuaU^ though we notice 
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straight away that we have scaling of the temperature 
even before we go to higher order terms. Thus, even 
though we are looking at a 2-D problem we have effec- 
tively higher dimensional behavior. As is knownia the 
Ising model is equivalent to a Coulomb gas in all dimen- 
sions and this is always in the plasma phase for d > 2. 
Our final results from the perturbative renormalization 
group to first order in g yield the following differential 
renormalization equations 



Finite Temperature RG flows 



ds 
d(m 2 ) 
ds 



r[ 



8m 2 y^ 2 / 3 



9, 



—m 

2 



(23) 
(24) 




The resultant flows in the {m 2 ,g} plane are shown in 

Fig. The differential RG equations show that m 2 al- FIG - 5: Renormalization group flows in {m 2 ,g} parameters. 



ways grows and that the flow of the variable g changes 
sign from positive to negative as m 2 becomes smaller 

than 207^573 ' indicated by point c on the plot. How- 
ever this point is of no special importance because there 
is no phase transition: recalling that m 2 ~ y we see 
that the temperature always flows to zero. Thus, the 
system remains in the same phase with flat domain walls 
that it has at T = 0, even as the temperature becomes 
finite. This formal RG result thus confirms the earlier 
expectation based on Eq. JTSJ. This result is to be con- 
trasted with the first order flow equations for the case 
of short range interactions where s^p- = and flows are 
vertical at this level of approximation - cf. Eq. (2.24) in 
Ref. ^3 truncated to first order in y, with y and J being 
the analogs of g and m 2 respectively in our analysis. It is 
interesting to note that simple tree level RG rescaling in 
a d-dimensional (short-range) sine-Gordon model yields 
4^ = (d — 2) J, in the analogous notation of Ref. [l8| . 
suggesting that our model is in effectively | dimensions. 

In addition to this classical case, we have also consid- 
ered the quantum mechanical case, which is equivalent 
to a classical sine-Gordon model in one extra dimension. 
Since this raises the dimensionality of the system even 
more, we find that there is no quantum roughening ei- 
ther. 



tigate the effect of the long range forces at finite tem- 
peratures. A similar problem has been investigated by 
Thouless 21 « who looked at one-dimensional Ising systems 
that have long range interactions ~ r~ 2 . He concluded 



that for interactions that fall off faster than 



there 



cannot be an ordered state, whereas for an interactions 
that falls off as r~ 2 the order cannot go continuously to 
zero. Even though our one-dimensional interface prob- 
lem is a height model, not an Ising model, the Thouless 
results suggest that the interface may be rough at finite 
temperatures. 

As before we obtain the energy for the step fluctuations 
in analogy with electrostatics, 



Us = 2m a / dzdz' 



h'{z)ti{z') 



(h(z) - h{z')f + (z - z') 2 



1/2- 



(26) 



By assuming that 6 < 1 we get to zeroth order 
in (h(z) — h{z')) the step fluctuation energy in Fourier 
space, 



■> / dk z dky 
Usq = rnm a ' 



( 2 ^) 2 (k 2 + k 2 ) 



B. Further Possibilities: Stepped Surfaces 

It is also interesting to consider the possibility of steps 
and investigate whether the steps themselves are smooth 
or rough. As our system we considered two semi-infinite 
three-dimensional domains of antiparallcl magnetization 
in the z-direction and added to the domain interface a 
unit lattice step in the x-direction. The magnetization 
vector in this case is defined by a one-dimensional inter- 
face height function h(z), 

M = {0, 0, m sgn {x - aQ[y - h(z)])} . (25) 

For the case of short-range forces as shown in Fig. 0] 
the classical interface is always rough. Here we inves- 



Our approximation scheme here introduces the necessity 
for an ultra-violet cutoff, A y in the k y integral, which is 
set by the discreteness of the lattice, and namely that the 
two height variables considered in the correlator are not 
evaluated at arbitrarily close points, but are at the very 
least on adjacent lattice sites. 

If we now consider the mean height correlation, and 
once more introduce a surface tension term, ~f s k 2 , we 
have, 



G s = ((h(z)-h(O)) 2 ) 



dk 7 . 



27r 4/3m 2 fc 2 ln(^ 



(28) 
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In the large distance limit G s ~ t!tl\ and the Gaus- 
sian model is rough at all finite temperatures. The long 
range interaction gives rise to a logarighmic term in the 
energy, raising the effective dimensionality of the system 
only very slightly above one. This suggests that if we 
were to take into account the lattice structure by adding a 
periodic pinning potential to form another modified sine- 
Gordon model there would be no effect on the roughness 
of the interface since the periodic sine term cannot pin 
the interface in dimensions lower than two. This is clear 
if we keep in mind that the usual sine-Gordon model for 
short range interactions only displays a phase transition 
in exactly two dimensions, where the free Gaussian model 
has logarithmically divergent correlations. In the usual 
sine-Gordon model in one dimension the sinusoidal pin- 
ning potential is unable to control the linearly diverging 
correlations of the free Gaussian model, and hence inter- 
faces remain rough. This argument suggests that for the 
long-range sine-Gordon model in one dimension where 
the effective dimensionality is close to one, the parame- 
ters flow to the high temperature limits, and just as for 
the case of short range forces, the periodic potential is 
unable to pin the step which fluctuates freely and hence 
is rough at all finite temperatures. 

For the quantum case with short range forces, G s di- 
verges as ~ In L for the Gaussian model and the full sine- 
Gordon model exhibits a roughening transition in the KT 
universality class as shown in Fig.0J Here, for the case of 
long-range forces in the quantum problem, similar anal- 
ysis shows that that the divergence of G s ~ Vln L is 
extremely weak for the Gaussian problem, possibly indi- 
cating that the sine-Gordon model will be smooth though 
we have not been able to verify this in a full RG calcula- 
tion. 
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APPENDIX: RENORMALIZATION DETAILS 

In the following we give some more details of the 
smooth cutoff procedural leading to the differential 
renormalization group equations l|23|l and l|24|l . Begin- 
ning with our choice for a smooth cutoff function, 



f(k y ,k z ) = e- a2 ^+^ 3 \ 



(A.l) 



and using k' z — e^ s , k' — e s k y , we obtain for small s a 
rescaled cutoff 



f'(k' y ,k' z ) = e ^ (A.2) 
« f(k' y ,k' z ) + ak' y ,k' z ) + 0( S 2 ), (A.3) 



with 



,[7 i ; / .A<i = :-W l (^ + 7 |A//)/(/,; / .A' : ). (A.4) 



We can now define a new field, say x(f), having £ as 
its smooth cutoff. We consider a sine-Gordon model for 
both the x an d </> fields and write 

S' = So W ,f} + Si [x, C] + Soi W, x] , ( A.5) 

with 



IV. DISCUSSION 



A; 



+ 



7T J {2n) 2 f(k y ,k z ) la\k y \ 
+l(k 2 y )}W{k)\\ (A.6) 



We have found that even though LiHoF4 interfaces ini- 
tially appear to be ideal for having a roughening transi- 
tion, the same long-range interactions which account for 
the system's domain structure turn out to be also re- 
sponsible for the lack of a roughening transition. The 
long-range interaction term which arises from the dipo- 
lar interactions effectively raises the dimensionality of the 
system and makes it equivalent to a sine-Gordon model 
in dimensions greater than two. The RG flow diagram, 
Fig. El which we obtain for m 2 ~ ^ indicates that in 
effect the roughening transition coincides with the bulk 
order-disorder transition and the whole ferromagnetic re- 
gion is smooth. In contrast, a step in the interface seems 
to have a rough profile at all temperatures, and the long- 
range interactions do not raise the dimensionality of the 
step enough to drive it into the smooth phase. 



m 2 f d 2 kez s r k 2 
liX ' CJ = ~~J (2w) 2 ak y ,k z )\^\k y ~\ + 

+7(^)}|x(fc)| 2 , (A.7) 

and 

SoiW,x] = d 2 rcos(0'(r) + X (r)). (A.8) 

It is straightforward to showi^ that the action S' re- 
sults in the same partition function corresponding to the 
rescaled action S in Eq. (TSJ. In Ref. [18| this is done 
by mapping to the Coulomb gas, but an easier method is 
simply to shift the argument of the cosine term by — xif) 
and then carry out the Gaussian integral over x( r )j to 
recover the rescaled action S. 



Now, as part of the RG calculation we integrate out 
the extra field \ to restore the original cutoff function. 
We use a cumulant expansion, in 5*oi and write 

/ D(j}'e So <e Sol ) 1 



where the averaging < >i indicates integrating out 
the extra x(r) fields using the Gaussian action S\. In 
addition, a multiplicative constant has been absorbed 
into the measure. 



The average of x( r ) 2 is given by 



(Xir) 2 ) 1 



d 2 k air\k y \e ^ s C ) (ky,k z ) 
(2tt) 2 2m 2 {k 2 + 1 a\k y \ 3 } 



4m 2 y / 7T7 2 / 3 



We gather our results to obtain 



(A.13) 



In this analysis we consider the first term in the expo- 
nential of the cumulant expansion 



{So,), = ( \d\ iL cos(^(r)+x(r)) 



d 2 r ^ cos(^'(r)) < cos(x(r)) >i 



■sm(4>'(r)) < sin(x(r)) >i 



(A.10) 



The average over the sine term is zero by symmetry. For 
the other term we can use the result for Gaussian inte- 
grals which states 



so we obtain 



Gaussian 



(A.ll) 



(SdH = I d 2 r(^) cos0'(r) < cos(x(r) >a 



d 2 r (j^J cos0'(r)e"^ x( " )2 )i. (A.12) 



(5, 



01/1 



d 2 r 



cos <j>'(r)e 8™2^f t 2/3 



(|§) ( d 2 r cos 4>'{r) + 0{s 2 ), (A.14) 



where we have introduced the renormalized coupling (/,,, 
which is related to the original <? coupling: 



(1 



21 
3 J 



r[ 



1 H2;^w^J s J 5 ; (A - 15) 

From the above relation we obtain the differential 
renormalization equation (j23(l . Note that m 2 receives 
no further corrections to this order in the cumulant ex- 
pansion and hence Eq. I|24|l is obtained directly from 
Eq. (29 . 
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